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PREFACE 


The  work  described  in  this  report  was  performed  under  Contract  No.  DAAA15-93-C-0070, 
Design  and  Development  of  Regenerable  Air  Filtration  Units  and  Integration  into  Host  Applications.  This 
work  was  started  in  September  1994  and  completed  in  May  1995. 

The  use  of  either  trade  or  manufacturers’  names  in  this  report  does  not  constitute  an  official 
endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes  of  advertisement. 

This  report  has  been  approved  for  public  release.  Registered  users  should  request  additional 
copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should  direct  such  requests  to 
the  National  Technical  Information  Service. 
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CATALYTIC  OXIDATION  OF  CYANOGEN  CHLORIDE  OVER 
A  MONOLITHIC  OXIDATION  CATALYST 


1.  INTRODUCTION 

Present  air  purification  systems  designed  for  removal  of  chemical  warfare  agents 
from  air  streams  are  based  solely  on  activated,  impregnated  carbon,  namely  ASC 
whetlerite.  While  these  filters  function  well  against  a  wide  range  of  chemical  agents,  they 
possess  several  shortcomings.  First,  the  carbon  filter  has  a  limited  capacity  for  agents 
which  are  removed  by  chemical  reaction  and  those  which  are  weakly  adsorbed.  Second, 
prolonged  environmental  exposure  has  been  shown  to  reduce  the  capacity  of  these  filters 
for  agents  which  are  removed  by  chemical  reaction1.  The  result  of  these  shortcomings  is 
to  impose  change-out  requirements  which  may  present  logistical  as  well  as  disposal 
burdens  to  the  user.  Catalytic  oxidation  is  an  alternative  air  purification  technology  which 
is  being  investigated  as  a  means  of  alleviating  the  above  mentioned  burdens. 

Catalytic  oxidation  has  been  proven  to  be  a  viable  means  of  controlling  industrial 
offgas  emissions2, J.  Offgas  streams  containing  chlorinated  organic  molecules  present  an 
additional  concern  in  that  the  chlorine  associated  with  these  compounds  has  the  potential 
to  poison/inhibit  the  catalyst2,4'9.  Studies  by  Rossin  et  al. 4,8,9  report  that  the  oxidation  of 
chlorine-containing  compounds  over  supported  platinum  catalysts  will  be  significantly 
inhibited  due  to  the  adsorption  of  product  HC1  onto  catalytic  sites.  The  observed 
inhibition  is;  however,  reversible.  Deactivation  resulting  from  the  oxidation  of  chlorine- 
containing  compounds  has  been  observed  by  several  investigators2,4'9.  For  example, 
Rossin  and  Farris8  report  a  decrease  in  the  conversion  of  1,500  ppm  chloroform  over  a  2% 
Pt/a-Al203  catalyst  from  92  to  83%  over  the  duration  of  a  65  hour  experiment.  Lester10; 
however,  claims  that  a  titania  based  catalyst  is  capable  of  oxidizing  chlorine-containing 
compounds  without  deactivation. 

The  catalytic  oxidation  of  nitrogen-containing  compounds  has  the  potential  to 
form  NOxh'1j.  Lester  and  Homeyer11  studied  the  destruction  of  a  number  of  amines  over 
a  monolithic  catalyst  which  presumably  contained  platinum.  In  all  cases,  the  NOx 
selectivity  was  significant,  especially  at  high  reaction  temperatures.  Rossin12,  in  a  study 
involving  the  oxidation  of  hydrogen  cyanide  over  a  platinum  monolith,  reported  that  while 
the  catalyst  could  readily  destroy  hydrogen  cyanide,  NOx  selectivities  were  unacceptable. 
Later,  Rossin  and  Campbell13  reported  a  catalytic  process  for  direct  oxidation  of  nitrogen- 
containing  compounds  without  the  formation  of  NOx.  The  catalyst  employed  by  the 
process,  whose  details  were  not  disclosed,  represents  a  significant  advantage  in  selectivity 
over  conventional  oxidation  catalysts. 

The  objective  of  this  study  was  to  develop  an  experimentally  based  reactor  design 
equation  which  would  be  capable  of  describing  the  oxidation  of  cyanogen  chloride  over  a 
wide  range  of  process  conditions.  The  approach  taken  in  meeting  this  objective  was  to 
first  record  isothermal  reactor  design  data  over  a  wide  range  of  temperatures, 
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concentrations  and  residence  times,  and  then  correlate  the  design  data  using  a  reactor 
design  equation  which  takes  into  account  all  necessary  transport  resistances  and  kinetic 
effects.  Also  included  in  this  effort  was  quantitatively  identifying  all  reaction  products 
formed  during  the  oxidation  of  cyanogen  chloride,  assessing  the  effects  of  water  on 
catalyst  performance,  and  assessing  the  significance  of  catalyst  deactivation. 


2.  EXPERIMENTAL  METHODS 

Materials'.  Cyanogen  chloride  was  obtained  from  Aberdeen  Proving  Ground. 
Analysis  revealed  the  compound  to  be  greater  than  99%  pure.  Neat  cyanogen  chloride 
was  metered  into  a  17  liter  pressure  vessel  and  diluted  with  either  air  or  helium  such  that 
the  concentration  of  cyanogen  chloride  in  the  vessel  was  1.0%  by  volume.  Calibration 
gases  (namely  CO,  C02,  NO  and  N02)  were  obtained  from  Matheson.  The  catalyst 
employed  during  this  study  was  the  Military  Air  Purification  (MAP)  catalyst  and  was 
obtained  from  Allied  Signal.  The  catalyst  was  a  platinum  monolith  with  a  cell  density  of 
93  cells/cm2. 

Catalyst  Preparation'.  Catalyst  cores  were  cut  from  a  30  cm  diameter  by  10.2  cm 
long  monolith  block  using  a  diamond  tip  hole  saw.  The  catalyst  was  a  ceramic  monolith 
possessing  93  channels  per  cm2.  Once  cut,  the  resulting  catalyst  core  was  approximately 
2.2  cm  diameter  by  6.0  cm  long.  The  center  channels  of  the  monolith  were  plugged  with 
alundum  cement  so  that  only  the  channels  near  the  outer  circumference  of  the  catalyst  core 
(totaling  115  channels)  remained  open.  The  resulting  catalyst  volume  was  7.42  cm3. 
Preparing  the  catalyst  in  this  manner  minimizes  axial  and  radial  temperature  gradients, 
thereby  allowing  one  to  record  reaction  rate  data  under  near  isothermal  conditions12,14. 
Following  preparation,  the  catalyst  core  was  wrapped  with  a  thin  layer  of  glass  wool  and 
loaded  into  the  reactor.  The  glass  wool  was  used  to  maintain  a  seal  between  the  catalyst 
and  the  reactor  wall,  thereby  preventing  the  gas  stream  from  by-passing  the  catalyst.  Once 
loaded  into  the  reactor,  seven  type  K  fine-wire  thermocouples  were  extended  axially  into 
unobstructed  channels  of  the  monolith.  These  thermocouples  were  located  0.5  cm  above 
the  monolith,  and  0.5,  1.0,  2.0,  3.0,  4.0  and  5.0  cm  from  the  inlet,  and  were  used  to  assess 
the  axial  temperature  gradient  within  the  monolith. 

Equipment:  A  schematic  representation  of  the  fixed  bed  reactor  system  is 
illustrated  in  Figure  1.  Dry,  oil-free  air  from  a  PSA  air  drier  was  metered  to  the  reactor 
using  either  a  0-5  or  0-20  Nl/min  mass  flow  controller.  N1  is  defined  as  one  liter  of  dry  air 
at  0°C  and  one  atmosphere  pressure.  Dry  air  from  the  mass  flow  controller  was  then 
delivered  to  a  water  saturator,  which  was  housed  in  a  temperature  controlled  chamber.  A 
back  pressure  regulator  is  located  downstream  of  the  water  saturator.  The  concentration 
of  water  in  the  air  stream  may  be  controlled  by  controlling  the  temperature  and  pressure  of 
the  water  saturator.  In  addition,  the  water  saturator  could  be  by-passed  to  conduct  the 
experiments  in  dry  air. 
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Figure  1:  Schematic  representation  of  Fixed  bed  catalytic  reactor  system.  BPR: 
back  pressure  regulator;  GC:  gas  chromatograph;  MFC:  mass  flow  controller. 


Cyanogen  chloride  was  delivered  to  the  air  stream  from  a  compressed  gas  cylinder 
containing  approximately  4.0%  cyanogen  chloride  in  helium.  The  reactant  stream  was 
metered  into  the  humidified  air  stream  via  one  of  three  mass  flow  controllers;  0-50 
Nml/min,  0-200  Nml/min  and  0-2,000  Nml/min.  Nml  refers  to  normal  milliliter  and  is 
defined  as  one  liter  of  air  at  0°C,  one  atmosphere  pressure.  In  some  cases,  the  0-2,000 
Nml/min  flow  controller  was  replaced  with  a  0-5,000  Nml/min  controller.  Use  of  this  flow 
controller  bank  allowed  for  achieving  a  wide  range  of  cyanogen  chloride  concentrations 
over  the  desired  air  flow  rate  range. 

Once  blended,  the  feed  stream  was  delivered  to  the  catalytic  reactor.  The  reactor 
consisted  of  a  66  cm  long  glass  tube  approximately  2.5  cm  in  diameter.  The  upper  40  cm 
of  the  glass  reactor  served  as  the  pre-heat  zone  and  was  filled  with  4  mm  diameter  glass 
beads.  The  glass  beads  provided  further  mixing  of  the  feed  gas  as  well  as  surface  area  for 
heating  the  incoming  feed  stream  to  reactor  temperature.  The  reactor  tube  was  housed  in 
an  8.9  cm  diameter  by  60  cm  long  aluminum  block.  The  aluminum  block  was  electrically 
heated,  and  the  temperature  of  the  reactor  was  controlled  by  controlling  the  temperature 
of  the  aluminum  block.  The  catalyst  core  was  located  approximately  7  to  10  cm  below  the 
pre-heat  zone  of  the  reactor. 

Following  the  reactor,  the  effluent  stream  was  delivered  to  a  catalytic  filter,  which 
was  designed  to  destroy  any  unreacted  cyanogen  chloride  in  the  effluent  stream. 
Following  the  catalytic  filter,  the  air  stream  was  delivered  to  two  scrubber  tanks  in  series. 
Each  scrubber  tank  was  filled  approximately  half  way  with  a  10  molar  sodium  hydroxide 
solution  and  served  to  remove  hydrochloric  acid  (a  reaction  product)  from  the  effluent 
stream.  Flow  from  the  scrubber  vessels  was  delivered  to  a  back  pressure  regulator,  which 
served  to  maintain  a  constant  pressure  on  the  reactor.  Following  the  back  pressure 
regulator,  the  air  stream  was  vented  into  a  fume  hood.  All  lines  in  contact  with  cyanogen 
chloride  were  contained  within  the  fume  hood  and  were  electrically  heat  traced. 

Isothermal  Reaction  Rate  Measurements.  Reaction  rate  data  were  recorded  at 
temperatures  of  200,  250,  275  and  300°C  at  a  pressure  of  6±1  psig.  The  feed 
concentration  of  cyanogen  chloride  was  varied  between  100  and  10,000  ppm  (274  and 
27,400  mg/m3)1  in  humid  (Tdew  =  23±1°C)  air.  The  residence  time,  based  on  the  reactor 
volume  occupied  by  the  unobstructed  monolith  channels  and  calculated  at  0°C  and  1  atm 
pressure,  was  adjusted  in  an  effort  to  achieve  conversions  between  10  and  95%.  Air  flow 
rates  were  typically  varied  between  2  and  25  Nl/min.  Reaction  rate  data  were  recorded 
only  during  day-time  hours.  Overnight,  the  catalyst  was  maintained  at  reaction 
temperature  under  flowing,  humid  air.  The  same  catalyst  core  was  used  for  all 
experiments  conducted  in  this  study.  All  process  conditions  were  maintained  for  between 
1  to  3  hours  to  ensure  the  achievement  of  steady-state.  It  should  be  noted  that  for  feed 
concentrations  greater  than  about  1,000  ppm,  an  increase  in  the  catalyst  temperature  was 
observed.  In  these  instances,  the  catalyst  temperature  was  adjusted  such  that  the  average 
catalyst  temperature  (numerical  average  based  on  the  thermocouple  measurements)  was  ± 
2°C  of  the  desired  value.  Both  feed  and  effluent  streams  were  analyzed  continuously  for 


1  Concentrations  reported  in  mg/m3  are  referenced  to  0°C,  one  atmosphere  pressure. 
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C02  and  cyanogen  chloride  during  all  runs.  The  conversion  of  cyanogen  chloride  was 
determined  based  on  the  concentration  of  cyanogen  chloride  in  the  feed  and  effluent 
streams. 

Effects  of  Water.  The  effects  of  water  concentration  on  the  catalytic  activity  were 
determined  by  recording  the  conversion  of  cyanogen  chloride  as  a  function  of  reaction 
temperature  for  water  concentrations  of  <0.03%,  0.1%,  0.25%,  0.60%,  1.14%  and 
2.71%.  These  concentrations  correspond  to  dew  point  temperatures  of  <-35°C,  -22°C, 
-11°C,  0°C,  10°C  and  23°C,  respectively,  at  one  atmosphere  pressure.  The  dew  point 
temperature  of  <-35°C  corresponds  to  air  delivered  from  the  PSA  air  drier.  The  run  was 
initiated  by  heating  the  reactor  to  320°C  in  flowing,  humid  air  at  a  residence  time  of  0.25 
seconds.  Once  at  temperature,  cyanogen  chloride  was  introduced  to  the  feed  stream  at  a 
rate  required  to  achieve  a  concentration  of  1,000  ppm.  These  process  conditions  were 
maintained  for  approximately  one  hour,  at  which  time,  the  catalyst  temperature  was 
decreased  to  170°C  at  30°C/hr,  with  the  effluent  stream  sampled  every  20  minutes  (10°C 
temperature  intervals)  for  C02,  N20  and  cyanogen  chloride  in  the  feed  and  effluent 
stream.  For  the  run  conducted  with  a  dew  point  temperature  of  <-35,  the  run  was 
initiated  by  heating  the  reactor  to  400°C. 

Catalyst  stability.  The  stability  of  the  catalyst  was  evaluated  by  exposing  it  to 
1,000  ppm  cyanogen  chloride  in  humid  (Tdew  =  23±C)  air  employing  a  residence  time  of 
0.05  seconds  at  275°C.  These  process  conditions  were  maintained  for  68  hours,  at  which 
time  the  run  was  terminated.  During  the  run,  the  effluent  stream  was  sampled  for  the 
concentration  of  C02,  N20  and  cyanogen  chloride  every  hour.  Just  prior  to  termination  of 
the  run,  the  effluent  stream  was  sampled  for  the  concentration  of  HC1  and  Cl2  using 
Drager  tubes. 

Upon  termination  of  the  run,  an  additional  test  was  performed  in  dry  air  (Tdew  < 
-35°C)  employing  a  residence  time  of  0.50  seconds  and  a  cyanogen  chloride  feed 
concentration  of  1,000  ppm.  This  run  was  conducted  at  400°C.  Prior  to  the  start  of  the 
run,  the  catalyst  was  exposed  to  the  dry  air  stream  overnight  in  order  to  remove  all  traces 
of  water  from  the  lines  as  well  as  equilibrate  the  catalyst  with  respect  to  the  concentration 
of  hydroxyls  on  the  surface.  The  run  was  then  initiated  the  next  morning  and  discontinued 
following  8  hours  of  operation.  Upon  termination  of  the  run,  the  catalyst  was  cooled  to 
275 °C,  and  the  previously  conducted  stability  run  was  performed. 

Sample  Analysis:  During  the  reaction  rate  studies,  the  concentration  of  cyanogen 
chloride  in  the  feed  and  effluent  stream,  and  the  concentration  of  C02  and  N20  in  the 
effluent  stream  were  monitored  using  a  gas  chromatograph.  All  analyses  were  performed 
using  a  Hewlett-Packard  5890  gas  chromatograph  (GC)  equipped  with  both  a  flame 
ionization  detector  (FID)  and  a  thermal  conductivity  detector  (TCD)  and  automatic 
sampling  valves.  A  sample  of  the  GC  trace  is  reported  in  Figure  2.  Cyanogen  chloride 
was  analyzed  using  a  2  m  by  3.2  mm  diameter  Hayesep  Q  column  attached  to  the  FID, 
while  C02  and  N20  were  analyzed  using  a  2  m  by  3.2  mm  diameter  Hayesep  Q  column 
attached  to  the  TCD.  Samples  were  analyzed  first  by  sampling  the  effluent  for  C02  and 
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Figure  2:  Sample  GC  trace  showing  feed  and  effluent  analysis.  1.  Air,  2.  CO2,  3.  N2O, 
4.  H2O,  5.  Cyanogen  Chloride  (effluent),  6.  Cyanogen  Chloride  (feed). 
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N2O  at  60°C,  then  heating  the  oven  to  120°C  at  40°C/min  and  sampling  the  effluent  for 
cyanogen  chloride.  Following  this,  the  feed  stream  was  sampled  for  cyanogen  chloride. 
The  above  analysis  sequence  required  17  minutes  to  perform. 


3.  MODEL  DEVELOPMENT  AND  FIT  PARAMETER  ESTIMATION 


Reaction  Rate  Model  Development.  The  monolithic  oxidation  catalyst  employed 
in  this  study  consisted  of  a  series  of  straight,  parallel  channels.  The  monolith  substrate  is 
coated  with  a  thin  layer  of  catalyst  support  material,  termed  washcoat.  The  catalytic 
metals  are  impregnated  onto  the  washcoat.  The  monolith  employed  in  this  study  had 
square  channels,  with  a  channel  density  93  cells/cm2.  The  thickness  of  the  washcoat  and 
substrate  were  not  supplied  by  the  manufacturer.  For  modeling  purposes,  the  washcoat 
was  assumed  to  be  0.0035  cm  thick,  and  the  substrate  was  assumed  to  be  0.01 15  cm  thick. 
A  schematic  representation  of  a  monolith  channel  is  provided  below  in  Figure  3.  Based  on 
these  dimensions,  the  monolith  is  1 1.55%  washcoat  plus  catalyst. 


.0922  crrrH 


Figure  3 :  Schematic  Representation  of  a  Monolith  Channel. 

The  oxidation  of  cyanogen  chloride  was  modeled  according  to  the  reaction  rate 
expression: 


Rate  = 


_ acck _ 

1  +  bchci  +  DCck 


(1) 


In  the  above  reaction  rate  expression,  A,  B  and  D  are  fit  parameters  (s'1,  cm3/mol  and 
cm3/mol,  respectively),  Chci  is  the  concentration  of  hydrochloric  acid  (mols/cm3),  Cck  is 
the  concentration  of  cyanogen  chloride  (mols/cm3)  and  Rate  is  the  reaction  rate  (mols/s- 
cm3  catalyst).  The  reaction  rate  expression  may  be  derived  by  assuming  the  reaction 
occurs  via  adsorption  of  cyanogen  chloride  onto  surface  hydroxyls  followed  by 
decomposition  of  the  adsorbed  species.  The  concentration  of  oxygen  was  held  constant 
during  all  experiments  (21%)  and  was  in  great  excess  (meaning  that  regardless  of  the 
cyanogen  chloride  conversion,  the  concentration  of  oxygen  did  not  change  significantly 
over  the  length  of  the  reactor).  Therefore,  the  concentration  of  oxygen  was  not  taken  into 
account.  The  rate  expression  also  does  not  take  into  account  the  concentration  of  water. 
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As  long  as  the  concentration  of  water  was  greater  than  the  concentration  of  cyanogen 
chloride,  the  concentration  of  water  did  not  affect  the  reaction  rate. 

The  above  rate  expression  was  incorporated  into  a  fixed  bed  reactor  design 
equation  which  took  into  account  external  mass  transfer  resistances  and  kinetic  effects. 
The  equations  which  govern  mass  transfer  and  chemical  reaction  rates  within  the  channel 
of  the  monolith  are15: 


d  C® 

v — ^  =  k 
dz 


m 


A[C 


s 

CK 


+ 


Cg  1 
'"CK-1 


(2) 


Rate  kmA  [Cg^  * 


(3) 


where  v  is  the  average  linear  velocity  within  the  channel  at  reaction  conditions  (cm/s), 

CgK  arid  CpK  refer  to  the  gas  and  solid  phase  concentration  of  cyanogen  chloride, 

respectively,  calculated  at  reaction  conditions  (mols/cm3),  dz  is  the  differential  reactor 
length  (cm),  km  is  the  mass  transfer  coefficient  (cm/s),  A  is  the  surface  area  of  the 
monolith  channel  (cm2/cm3  void),  Rate  is  the  reaction  rate  as  defined  by  eq.  1  (mols/s- 
cm3-cat),  and  A1  is  the  surface  area  of  the  catalyst  (cm2/cm3  catalyst).  The  above  reaction 
rate  expressions  were  numerically  integrated  using  a  multistep  method.  Data  were 
correlated  over  the  entire  data  set  simultaneously,  rather  than  by  the  standard  technique  of 
determining  the  fit  parameters  at  each  reaction  temperature,  followed  by  correlating  the 
data  using  Arrhenius  type  equations.  Therefore,  both  the  activation  energy/heats  of 
adsorption  and  pre-exponential  factors  were  determined  simultaneously.  The  algorithm 
used  to  determine  the  fit  parameters  employed  numerical  derivatives  and  minimized  the 
error  between  the  predicted  and  experimental  conversions.  The  mass  transfer  coefficient 
was  determined  from  a  limiting  Sherwood  number  correlation15'16: 


Sh  =  3.66(1  +  0.95ReSc-)0'45 

z 


(4) 


k-m 


2Rh 


(5) 


where  Sh  is  the  Sherwood  number,  Re  is  the  Reynolds  number,  Sc  is  the  Schmidt  number, 
d  is  channel  diameter  (cm),  Da^  is  the  diffusion  coefficient  for  cyanogen  chloride  in  air 
(cm2/s),  and  R^  is  the  hydraulic  radius  of  the  channel  (cm).  The  diffusion  coefficient  was 
estimated  to  be  0.067  cm2/s  at  0°C  and  was  allowed  to  increase  proportional  to 
temperature  to  the  power  of  1.5  (T15).  Using  the  reactor  design  equations,  the  fit 
parameters  were  determined  to  be: 

A  =  1.8631(08)  exp  (-6,075/T)  s'1 
B  =  1.9446(04)  exp  (+4,456/T)  cm3/mol 
D  =  3.0643(05)  exp  (+2,695/T)  cm3/mol 
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4.  RESULTS  AND  DISCUSSION 


Axial  Temperature  Deviation :  Oxidation  reactions  are  highly  exothermic.  When 
recording  reaction  rate  data,  it  is  desired  to  minimize  heat  transfer  effects  so  that  near 
isothermal  operation  of  the  catalyst  is  achieved.  In  this  manner,  kinetic  effects  will  not  be 
masked  by  heat  transfer  effects.  For  pelleted  and  granular  catalysts,  temperature 
excursions  in  the  laboratory  reactor  can  be  and  often  are  controlled  by  diluting  the  catalyst 
bed  with  inert  particles  of  a  similar  mesh  size.  The  monolithic  oxidation  catalyst  consists 
of  a  single  piece  of  substrate  and  thus  cannot  be  diluted  in  the  same  manner  as  granular 
catalysts  are.  Therefore,  recording  reaction  rates  under  near  isothermal  conditions 
presents  a  challenge  to  the  investigator. 

A  differential  mode  of  operation  presents  one  option  for  minimizing  the  axial 
temperature  gradients.  Stock  and  Lowe17,18  have  measured  quality  reaction  rate  data  for 
the  oxidation  of  carbon  monoxide  by  operating  the  monolith  under  differential  conditions. 
However,  a  differential  mode  of  operation  was  felt  to  be  inappropriate  for  the  present 
study.  This  is  because  the  low  residence  times  required  to  maintain  differential  conditions 
may  result  in  entrance  and  exit  effects  dominating  the  fluid  dynamics  within  the  monolith 
channel,  rather  than  the  laminar  flow  that  would  be  encountered  by  the  application.  Also, 
reaction  rate  inhibition  due  to  the  adsorption  of  reaction  products  onto  catalytically 
reactive  sites,  as  has  been  reported  previously  for  the  oxidation  of  chlorine-containing 
compounds8’9,  may  be  difficult  to  assess,  since  the  concentration  of  the  reaction  products 
will  be  minimal.  Reaction  products  could  always  be  added  to  the  feed  stream,  but  this 
would  only  add  an  additional  order  of  complexity  to  an  already  complex  problem. 

The  oxidation  of  cyanogen  chloride  in  air  is  highly  exothermic.  The  heat  of 
reaction  (at  298  K)  was  estimated  to  be  107.9  kcals/mol.  The  adiabatic  temperature  rise 
(corresponding  to  100%  destruction  of  CK)  for  cyanogen  chloride  concentrations  of 
1,000,  3,000  and  10,000  ppm  were  calculated  to  be  13.8,  41.4  and  138°C,  respectively. 
Because  of  the  heat  generated  by  the  reaction,  one  can  anticipate  difficulties  in  maintaining 
a  uniform  temperature  along  the  catalyst  length,  especially  as  the  concentration  of  reactant 
and  reaction  temperature  are  increased.  The  technique  to  minimize  the  axial  temperature 
gradient,  as  described  in  the  Experimental  Methods  section,  involved  limiting  the  heat 
generated  within  the  volume  of  reactor  occupied  by  the  monolith  core  by  obstructing  the 
center  channels  of  the  monolith.  Therefore,  the  reaction  would  only  proceed  in  channels 
located  near  the  reactor  wall,  and  the  ratio  of  the  heat  generated  to  surface  area  of  the 
reactor  tube  is  maximized. 

Figure  4  reports  the  temperature  deviation  for  the  oxidation  of  cyanogen  chloride 
at  300°C.  Temperature  deviation  refers  to  the  difference  between  the  catalyst  temperature 
recorded  at  discrete  points  within  the  monolith  and  the  average  of  the  temperature 
measurements.  Data  correspond  to  cyanogen  chloride  feed  concentrations  of  1,000,  3,000 
and  10,000  ppm  (2,740,  8,220  and  27,400  mg/m3)  and  were  recorded  under  conditions 
where  the  conversion  of  cyanogen  chloride  was  greater  than  90%.  As  expected,  the 
temperature  uniformity  of  the  catalyst  was  greatly  improved  at  lower  conversions  and  at 
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Dimensionless  Catalyst  Length 

Figure  4:  Axial  temperature  deviation  recorded  during  the  oxidation  of  cyanogen  chloride  at 
300°C. 

lower  reaction  temperatures.  For  feed  concentrations  of  3,000  ppm  (8,220  mg/m3)  and 
below,  the  axial  temperature  deviation  was  typically  less  than  ±5°C.  For  feed 
concentrations  of  10,000  ppm  (27,400  mg/m3),  the  axial  temperature  deviation  increased 
to  almost  ±10°C.  These  results  attest  to  the  near  isothermal  operation  achieved  via 
modification  of  the  monolith  core. 

Effects  of  Water-.  Figure  5  reports  conversion  as  a  function  of  temperature  for  the 
oxidation  of  1,000  ppm  (2,740  mg/m3)  cyanogen  chloride  in  air  with  different 
concentrations  of  water.  Data  reported  in  this  figure  were  recorded  at  a  residence  time  of 
0.25  s.  Results  presented  in  this  figure  show  that  increasing  the  concentration  of  water 
from  0  to  0.25%  (Tdew  =  <-35  and  -11°C)  has  a  significant  effect  on  the  conversion  of 
cyanogen  chloride;  however,  further  increasing  the  concentration  of  water  to  2.7%  (Tdew  = 
23 °C)  did  not  affect  the  conversion  of  cyanogen  chloride.  These  results  can  be  explained 
assuming  the  reaction  is  initiated  by  the  hydrolysis  of  cyanogen  chloride,  rather  than  by  an 
oxidative  mechanism.  A  hydrolysis  reaction  involving  cyanogen  chloride  would  proceed 
via  an  interaction  between  a  cyanogen  chloride  molecule  and  a  hydroxyl  group  on  the 
surface  of  the  catalyst  to  form  (possibly)  an  adsorbed  -OCN  species  plus  hydrochloric 
acid.  The  -OCN  species  may  then  decompose  to  form  the  reaction  products. 


16 


110 


100 

90 


80 


5$  70 
O  60 

CO 


I  50 
c 


30 


20 


10 


0  L 

0 


♦#***#XXXXXX  x 

♦  A  _ 


■ 

[H201  =  0.00  % 

** 

• 

[H201  =  0.10  % 

• 

• 

■V  ♦ 

[CKl 

=  1,000  ppm 

A 

[H20I  =  0.25  % 

• 

+  • 

tau 

=  0.25  s 

♦ 

[H201  =  0.60  % 

♦ 

►X 

• 

+ 

[H201  =  1.14  % 

• 

+ 

X 

IH201  =  2.71  % 

♦  x 

c  onv.pl  o 


50  100  150  200  250  300  350  400  450 

Temperature,  C 


Figure  5:  Conversion  as  a  function  of  temperature  for  the  oxidation  of  1,000  ppm  cyanogen 
chloride  in  air  as  a  function  of  water  concentration. 

One  possible  explanation  for  the  effects  of  water  on  the  catalyst  performance  stems 
from  the  chemisorption  of  water  on  the  surface  of  the  catalyst.  At  low  water 
concentrations,  the  surface  of  the  catalyst  is  not  saturated  with  hydroxyls;  i.e.  there  are 
sites  available  for  adsorption  of  water.  Under  these  conditions,  increasing  the 
concentration  of  water  in  the  feed  stream  will  increase  the  hydroxyl  concentration  at  the 
surface  of  the  catalyst,  thereby  increasing  the  catalytic  activity.  Upon  further  increasing 
the  concentration  of  water  in  the  feed  stream,  the  catalyst  surface  becomes  saturated  with 
hydroxyls;  i.e.  there  are  no  more  sites  available  for  water  adsorption.  As  a  result,  the 
reaction  rate  is  no  longer  a  function  of  the  concentration  of  water  in  the  feed  stream. 

Reaction  Products :  The  only  carbon-containing  reaction  product  identified  in  the 
reactor  effluent  stream  was  carbon  dioxide.  No  carbon  monoxide  or  products  of  partial 
oxidation  were  identified.  Carbon  balances  were  typically  100±5%  for  all  runs  conducted 
in  this  study.  Figure  6  reports  the  product  distribution  of  the  nitrogen-containing  reaction 
products  as  a  function  of  temperature  for  the  oxidation  of  1,000  ppm  cyanogen  chloride  at 
atmospheric  pressure.  This  test  was  conducted  by  placing  a  small  piece  of  monolith  into  a 
0.95  cm  o.d.  glass  tube  and  performing  the  test  in  the  micro-reactor  system.  In  this 
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Figure  6:  Nitrogen-product  concentration  (A)  and  nitrogen-product  selectivity  (B)  as  a  function  of 
reaction  temperature  for  the  oxidation  of  1,000  ppm  cyanogen  chloride  in  humid  air. 


manner,  near  isothermal  operation  could  be  achieved  and  the  product  distribution  would 
not  be  skewed  as  a  result  of  temperature  effects.  Nitrogen-containing  reaction  products 
consisted  of  N2,  N2O,  NO  and  NO2,  with  NO  plus  NO2  being  collectively  referred  to  as 
NOx  in  Figure  6.  The  product  distribution  was  found  to  vary  as  a  function  of  temperature, 
with  the  formation  of  NOx  favored  at  high  temperatures,  and  the  formation  of  the  reduced 
products  (N2  and  N20)  favored  at  the  lower  temperatures.  At  400°C,  approximately  65% 
of  the  nitrogen-containing  products  (on  an  atomic  nitrogen  basis)  consist  of  NOx.  Below 
300°C,  NOx  decreases  to  less  than  15%;  however,  under  no  conditions  where  the 
cyanogen  chloride  conversion  was  greater  than  99%  was  NOx  not  measured  in  the  effluent 
stream.  Results  presented  in  Figure  6  are  important  because  NOx  will  be  difficult  to 
remove  from  a  reactor  effluent  stream.  Campbell  and  Rossin  (1995)  reported  a  catalyst 
capable  of  destroying  nitrogen-containing  compounds  without  the  formation  of  NOx. 
Such  a  catalyst  would  be  required  to  meet  the  demands  associated  military  chemical 
defense  applications. 


Catalytic  Deactivation.  Figure  7  reports  the  conversion  of  1,000  ppm  (2,740 
mg/m3)  cyanogen  chloride  as  a  function  of  time-on-stream  in  humid  air.  Data  were 
recorded  at  275°C  employing  a  residence  time  of  0.05  seconds.  Results  illustrate  that  the 


Figure  7:  Conversion  as  a  function  of  time-on-stream  for  the  oxidation  of  1,000  ppm  cyanogen 
chloride  in  humid  air  at  275°C. 
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conversion  of  cyanogen  chloride  remain  constant  at  about  92%  throughout  the  duration  of 
the  run.  During  the  run,  7.4  cm3  of  monolith  (about  0.85  g  of  washcoat)  destroyed  99.6  g 
of  cyanogen  chloride.  Results  demonstrate  that  the  catalyst  is  stable  and  will  not  be 
deactivated  by  exposure  to  cyanogen  chloride. 

Following  completion  of  the  run,  the  catalyst  was  exposed  to  dry  air  (Tdew  =  <- 
35°C)  at  400°C  overnight  in  order  to  remove  all  traces  of  moisture  from  the  system.  The 
stability  of  the  catalyst  was  then  evaluated  under  dry  conditions  by  exposing  the  catalyst  to 
1,000  ppm  cyanogen  chloride  at  400°C  at  a  residence  time  of  0.5  s.  The  results  of  this 
test  are  shown  in  Figure  8.  The  conversion  of  cyanogen  chloride  rapidly  decreased  over 
the  first  hour  of  the  run.  Following  this  time,  the  conversion  continued  to  decrease,  but  at 
a  slower  rate.  A  comparison  of  results  presented  in  Figures  7  and  8  highlights  the ’effects 
of  water  on  the  performance  of  the  catalyst,  In  dry  air,  increasing  the  reaction 
temperature  by  125°C  and  decreasing  the  residence  time  by  a  factor  often  was  required  in 
order  to  achieve  a  meaningful  conversion.  Following  this  run,  the  catalyst  was  cooled  to 
275°C  overnight  in  flowing,  humid  air.  In  the  morning,  the  catalyst  was  exposed  to  1,000 
ppm  cyanogen  chloride  at  275°C  at  a  residence  time  of  0.05  s  in  humid  air.  The 
conversion  of  cyanogen  chloride  was  approximate  90%,  consistent  with  the  conversion 
obtained  prior  to  the  dry  air  exposure.  Interestingly,  exposure  of  the  catalyst  to  cyanogen 
chloride  under  dry  conditions  did  not  irreversibly  deactivate  the  catalyst. 


Time,  hrs 

Figure  8:  Conversion  as  a  function  of  time-on-stream  for  the  oxidation  of  1,000  ppm  cyanogen 
chloride  in  dry  air  at  400°C. 
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Reaction  Rate  Measurements'.  Axial  temperature  gradients  were  assessed  during 
each  run  from  temperature  measurements  obtained  using  thermocouples  located  within 
channels  near  the  outer  circumference  of  the  monolith.  Typically,  the  axial  temperature 
deviation  was  less  than  3-5°C,  with  the  deviation  increasing  to  10°C  for  runs  conducted 
with  a  10,000  ppm  feed  concentration.  Even  the  10°C  temperature  deviation  was  not 
deemed  to  be  unacceptable,  and  therefore  data  recorded  under  these  conditions  were  used 
in  the  fit  parameter  analysis. 

Reaction  rate  data  were  recorded  at  temperatures  of  200,  250,  275  and  300°C  for 
cyanogen  chloride  feed  concentrations  between  100  and  10,000  ppm  (274  and  27,400 
mg/m3)  at  6±1  psig  in  humid  (Tdew  =  23±1°C)  air.  Figures  9  through  12  report  the 
conversion  of  cyanogen  chloride  as  a  function  of  residence  time  at  reaction  temperatures 
of  200,  250,  275  and  300°C,  respectively.  Results  presented  in  these  figures  show  that 
increasing  the  concentration  of  cyanogen  chloride  has  a  significant  effect  on  conversion, 
especially  at  200°C.  For  example,  increasing  the  feed  concentration  of  cyanogen  chloride 
from  300  ppm  to  1,000  ppm  at  200°C  at  a  residence  time  of  0.1  seconds  results  in  the 
conversion  decreasing  from  about  75%  to  about  50%.  At  higher  reaction  temperatures, 
the  effects  of  increasing  the  cyanogen  chloride  feed  concentration  on  conversion  are 
present  but  less  pronounced.  Results  indicate  that  the  reaction  is  being  inhibited  by  the 
adsorption  of  reactant  and/or  reaction  product  (HC1)  on  catalytic  sites. 

Several  different  forms  of  the  reaction  rate  expression  were  evaluated  for  their 
ability  to  correlate  the  experimental  data.  The  form  of  the  reaction  rate  expression  which 
best  described  the  experimental  data  is  reported  as  eq.  1 .  This  rate  expression  takes  into 
account  the  adsorption  of  cyanogen  chloride  and  product  hydrochloric  acid  on  the 
reaction  sites.  Failure  to  account  for  the  adsorption  of  hydrochloric  acid  resulted  in  a  poor 
data  correlation.  The  solid  lines  in  the  Figures  9  through  12  represent  the  data  correlated 
using  the  reactor  design  equation  (eqs.  2  and  3).  Figure  13  shows  a  parity  plot  of  the 
predicted  and  experimental  conversion  of  cyanogen  chloride.  Results  presented  in  Figures 
9  through  13  demonstrates  the  ability  of  the  design  equation  to  accurately  describe  the 
experimental  data. 

The  form  of  the  reaction  rate  expression  and  the  effects  of  water  concentration 
(Figure  5)  suggest  that  the  reaction  is  occurring  via  a  catalyzed  hydrolysis  reaction,  with 
the  reaction  being  inhibited  by  the  presence  of  product  hydrochloric  acid.  The  reaction 
rate  expression  may  be  derived  by  assuming  gas  phase  cyanogen  chloride  adsorbs  onto 
surface  hydroxyls  and  is  subsequently  decomposed,  with  the  reaction  being  inhibited  by 
the  presence  of  HC1.  It  should  be  noted  that  the  reaction  rate  expression  did  not  take  into 
account  the  effects  of  water.  This  is  because  data  used  to  derive  the  rate  expression  were 
recorded  under  conditions  in  which  the  concentration  of  water  did  not  affect  the  catalytic 
activity. 
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Figure  9:  Conversion  of  cyanogen  chloride  as  a  function  of  residence  time  in  humid  air  at 
-  200°C.  Solid  lines  represent  data  as  correlated  using  the  reactor  design 

equation. 


Figure  10:  Conversion  of  cyanogen  chloride  as  a  function  of  residence  time  in  humid  air 
at  250°C.  Solid  lines  represent  data  as  correlated  using  the  reactor  design 


Figure  1 1 :  Conversion  of  cyanogen  chloride  as  a  function  of  residence  time  in  humid  air 
at  275°C.  Solid  lines  represent  data  as  correlated  using  the  reactor  design 
equation. 


100 


80 
I  60 

L. 

> 

o  40 
O 

20 


0 

.001  .003  .01  .03  .1  .3  1 

Residence  Time,  s 


Figure  12:  Conversion  of  cyanogen  chloride  as  a  function  of  residence  time  in  humid  air 
at  300°C.  Solid  lines  represent  data  as  correlated  using  the  reactor  design 
equation. 
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Figure  13:  Parity  plot  of  predicted  versus  experimental  conversion  of  cyanogen  chloride. 

Figure  14  reports  the  conversion  of  cyanogen  chloride  as  a  function  of  residence 
time  at  300°C.  Data  were  recorded  for  feed  concentrations  of  3,000  and  10,000  ppm 
(9,133  and  27,400  mg/m3)  and  are  reported  for  the  purpose  of  evaluating  the  ability  the 
model  to  predict  conversions  greater  than  99%.  The  high  conversion  scenario  was 
evaluated  because  military  air  purification  systems  require  that  the  challenge  be  reduced  by 
several  orders  of  magnitude.  The  solid  lines  in  Figure  14  represent  the  cyanogen  chloride 
conversion  predicted  using  the  reactor  design  equation.  Results  presented  in  this  figure 
demonstrate  that  the  reactor  design  equation  can  be  used  to  accurately  predict  cyanogen 
chloride  reduction  ratios  up  to  four  orders  of  magnitude  (which  was  the  detection  limits  of 
our  instruments)  under  isothermal  conditions. 


5.  CONCLUSIONS 

Cyanogen  chloride  may  be  readily  oxidized  over  the  monolithic  oxidation  catalyst 
at  temperatures  greater  than  about  250°C.  Reaction  products  consisted  of  C02,  HCI,  N2, 
N20  and  NOx,  with  the  nitrogen  product  selectivity  being  a  strong  function  of  catalyst 
temperature.  Catalytic  deactivation  was  not  observed  under  humid  conditions,  indicating 
that  the  catalyst  is  very  stable.  In  the  absence  of  water,  catalyst  deactivation  was 
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Figure  14:  Reduction  ratio,  C  /  C®,  of  cyanogen  chloride  as  a  function  of  residence  time  at 

a  d 

300°C. 

observed.  For  water  vapor  concentrations  less  than  0.10%,  the  presence  of  water  in  the 
feed  had  a  significant  effect  on  the  activity  of  the  catalyst.  In  the  absence  of  water,  the 
oxidation  of  cyanogen  chloride  was  very  slow,  requiring  reaction  temperatures  greater 
than  400°C  in  order  to  proceed.  The  addition  of  a  small  amount  of  water  to  the  feed 
stream  greatly  enhanced  the  catalytic  activity.  For  water  vapor  concentrations  greater 
than  0.25%,  the  concentration  of  water  in  the  feed  stream  did  not,  to  within  experimental 
error,  affect  the  catalytic  activity.  Under  isothermal  conditions,  the  oxidation  of  cyanogen 
chloride  could  be  modeled  employing  a  reaction  mechanism  which  assumes  the  reaction 
occurs  via  adsorption  and  decomposition  of  cyanogen  chloride  onto  surface  hydroxyl 
groups.  The  reactor  design  equation  was  able  to  accurately  describe  the  conversion  of 
cyanogen  chloride  over  a  wide  range  of  process  conditions. 
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Appendix 

FORTRAN  routine  used  to  determine  kinetic  fit  parameters 


c - 

C  This  routine,  MONOl.F,  determines  the  kinetic  fit  parameters 
C  for  the  oxidation  of  compound  in  air  over  a  monolith 
C  with  a  reaction  rate  expression: 

C 

C  ACa 

C  Rate  = - 

C  (1+BCp  +  DCa) 

C  Where: 

C  A,  B  and  D  are  fit  parameters 
C  Ca  is  the  reactant  concentration 
C  Cp  is  the  concentration  of  reaction  product 
C 

C  The  reaction  is  modeled  according  to: 

C 

C  v  dCa/dz  =  kA(Ca,g  -  Ca,s) 

C  Rate  =  kA(Ca,s  -  Ca,g) 

C 

C  Where  v  is  the  velocity,  z  is  the  reactor  length,  k  is 
C  the  mass  transfer  coefficient,  and  A  is  the  surface  area 
C  of  the  monolith.  The  subscripts  s  and  g  refer  to  the 
C  gas  and  solid  phase. 

C - 

IMPLICIT  REAL*8  (A-H.O-Z) 

DIMENSION  Y0(5),  CONV(IOO),  Al(6,7),  P5(6) 

DIMENSION  CONX(30),  X(100) 

COMMON  /SBCOM1/  XKM1,  SA,  SS,  CAO,  RHOB 
COMMON  /SBCOM3/  XI,  X2,  X3 
COMMON  /SBCOM2/  WASHT,  CELLR,  TOUT,  VEL 
common  /sbcom3/  j,  j5 

DATA  ONE  /1.0D00/,  TWO  /2.0D00/,  FOUR  /4.ODO0/,  ZERO  /0.0D00/ 
DATA  RG  /82.085D00/,  KORD  111,  KORD1 131 


N  is  the  number  of  equations  to  be  solved  by  EPISODE 


N  =  1 

OPEN(9  ,FILE='stat.  out') 

OPEN(7,FILE='test.out') 

OPEN(5,  FILE- mono  1  .dat') 

REWIND(5) 

C - 

C  NDAT  is  the  number  of  data  points. 

C  AO,  BO  and  DO  are  initial  guesses  for  the  pre-exponential  factors 
C  corresponding  to  the  fit  parameters  XI,  X2  and  X3.  EA,  EB  and 
C  ED  are  initial  guesses  for  the  corresponding  activation  energies. 

C - - - 

READ(5,28)  NDAT,  AO,  BO,  DO,  EA,  EB,  ED 
28  FORMAT(I7,/,D18.8,/,D18.8,/,D18.8,/,D18.8,/,D18.8,/,D18.8) 
READ(5,26)  RHOB,  SAVEH,  TOUT,  EPS 
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r>  n  n  o 


READ(5,27)  TREF,PRESS, CELLR,  WASHT,SUBT,XKM 
RE  AD(5,27)  DELA, DELB, DELD, DEL  1  ,DEL2,DEL4 

26  FORMAT(3(D18.8,/),D18.8) 

27  FORMAT(5(D18.8,/),D18.8) 

CL0SE(5) 

DELA  =  AO/DELA 
DELB  =  BO/DELB 
DELD  =  DO/DELD 
DELI  =E  A/DEL  1 
DEL2  =  EB/DEL2 
DEL4  =  ED/DEL4 


SS  surface  area,  cm**2  catalyst/cm**3  reactor  void 
SA  surface  area,  cm* *2  catalyst/cm**3  catalyst 


CELLR  =  CELLR  -  SUBT 
CELLT  =  CELLR  -  TWO*WASHT 
SS  =  FOUR/CELLT 

SA  =  FOUR*CELLT/(CELLR*CELLR  -  CELLT* CELLT) 
VOID  =  (CELLT*CELLT)/((CELLR+SUBT)*(CELLR+SUBT)) 

c - 

c  The  DO  600  loop  is  used  to  re-evaluate  the  function  with  a 
c  new  set  of  fit  parameters 

c - 

DO  600  L=  1,50 
FI  = ZERO 
F2  =  ZERO 
F3  =  ZERO 
F4  =  ZERO 
F5  =  ZERO 
F6  =  ZERO 
FI  A  =  ZERO 
FIB  =  ZERO 
FID  =  ZERO 
FI  1  =  ZERO 
F12  =  ZERO 
F14  =  ZERO 
F2A  =  ZERO 
F2B =  ZERO 
F2D  =  ZERO 
F2 1  =  ZERO 
F22  =  ZERO 
F24  =  ZERO 
F3A  =  ZERO 
F3B  =  ZERO 
F3D  =  ZERO 
F3 1  =  ZERO 
F32  =  ZERO 
F34  =  ZERO 
F4A  =  ZERO 
F4B = ZERO 
F4D =  ZERO 
F4 1  =  ZERO 
F42  =  ZERO 
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F44  =  ZERO 

F5A  =  ZERO 

FSB  =  ZERO 

F5D  =  ZERO 

F5 1  =  ZERO 

F52  =  ZERO 

F54  =  ZERO 

F6A  =  ZERO 

F6B  =  ZERO 

F6D  =  ZERO 

F6 1  =  ZERO 

F62  =  ZERO 

F64  =  ZERO 

OPEN  (2  ,FELE='ck.  dat') 

revvind(2) 

c - 

c  The  DO  500  loop  calculates  conversions  at  a  fixed  set  of  process 
c  conditions,  varying  the  fit  parameters  accordingly. 

c - 

DO  500  J  =  1.NDAT 

READ(2,12)  TAU,  TEMP,  CONC,  CONV(J) 

XKM1  =  XKM*(TEMP/573.15D00)**1.5D00 

VEL  =  TOOT*TEMP/(VOED*273. 15D00*PRESS*TAU) 

CAO  =  CONC 

TX  =  (ONE/TEMP  -  ONE/TREF) 

12  FORMAT(4D12.6) 

DO  400  J5  =  1,28 
IF(J5.EQ.l)  THEN 
XI  =  A0*DEXP(-EA*TX) 

X2  =  B0*DEXP(EB*TX) 

X3  =  D0*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ,2)  THEN 
XI  =  (A0  +  DELA)*DEXP(-EA*TX) 

X2  =  B0*DEXP(EB*TX) 

X3  =  D0*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.3)  THEN 
XI  =  (A0  -  DELA)*DEXP(-EA*TX) 

X2  =  B0*DEXP(EB*TX) 

X3  =  D0*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.4)  THEN 
XI  =  A0*DEXP(-EA*TX) 

X2  =  (B0  +  DELB)*DEXP(EB*TX) 

X3  =  D0*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.5)  THEN 
XI  =  A0*DEXP(-EA*TX) 

X2  =  (B0  -  DELB)*DEXP(EB*TX) 
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X3  =  DO*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ.6)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  (DO  +  DELD)*DEXP(ED*TX) 
ELSE 
END  IF 

IF(J5.EQ.7)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  (DO  -  DELD)*DEXP(ED*TX) 
ELSE 
END  IF 

IF(J5.EQ.8)  THEN 
XI  =  AO  *DEXP(-(E  A+DEL 1 )  *TX) 
X2  =  BO  *DEXP(EB  *TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ.9)  THEN 
XI  =  AO  *DEXP(-(E  A-DEL 1  )*TX) 
X2  =  BO*DEXP(EB*TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ.10)  THEN 
XI  =  AO  *DEXP(-E  A*TX) 

X2  =  B0*DEXP((EB+DEL2)*TX) 
X3  =  DO  *DEXP(ED  *TX) 

ELSE 
END  IF 

IF(J5.EQ.  1 1)  THEN 
XI  =  AO  *DEXP(-E  A*TX) 

X2  =  B0*DEXP((EB-DEL2)*TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.  12)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  D0*DEXP((ED+DEL4)*TX) 
ELSE 
ENDIF 

IF(J5.EQ.  13)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  BO  *DEXP(EB  *TX) 

X3  =  D0*DEXP((ED-DEL4)*TX) 
ELSE 
ENDIF 

IF(J5.EQ.  14)  THEN 
XI  =  (AO  +  DELA)*DEXP(-EA*TX) 
X2  =  (BO  +  DELB)*DEXP(EB*TX) 
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X3  =  DO*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.  15)  THEN 
XI  =  (AO  +  DELA)*DEXP(-EA*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  (DO  +  DELD)*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.  16)  THEN 

XI  =  (AO  +  DELA)*DEXP(-(E  A+DEL 1  )*TX) 
X2  =  BO*DEXP(EB*TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.  17)  THEN 
XI  =  (AO  +  DELA)*DEXP(-EA*TX) 

X2  =  B0*DEXP((EB+DEL2)*TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.18)  THEN 
XI  =  (AO  +  DELA)*DEXP(-EA*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  D0*DEXP((ED+DEL4)*TX) 

ELSE 

ENDIF 

IF(J5.EQ.  19)  THEN 
XI  =  AO  *DEXP(-E  A*TX) 

X2  =  (BO  +  DELB)*DEXP(EB*TX) 

X3  =  (DO  +  DELD)*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.20)  THEN 
XI  =  AO  *DEXP(-(E  A+DEL 1  )*TX) 

X2  =  (BO  +  DELB)*DEXP(EB*TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5.EQ.21)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  (BO  +  DELB)*DEXP((EB+DEL2)*TX) 
X3  =  DO*DEXP(ED*TX) 

ELSE 

ENDIF 

IF(J5EQ.22)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  (BO  +  DELB)*DEXP(EB*TX) 

X3  =  D0*DEXP((ED+DEL4)*TX) 

ELSE 

ENDIF 

IF(J5.EQ.23)  THEN 
XI  =  AO*DEXP(-(EA+DELl)*TX) 

X2  =  BO*DEXP(EB*TX) 
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X3  =  (DO  +  DELD)*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ.24)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  B0*DEXP((EB+DEL2)*TX) 

X3  =  (DO  +  DELD)*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ.25)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  (DO  +  DELD)*DEXP((ED+DEL4)*TX) 
ELSE 
END  IF 

IF(J5.EQ.26)  THEN 
XI  =  AO*DEXP(-(EA+DELl)*TX) 

X2  =  BO*DEXP((EB+DEL2)*TX) 

X3  =  DO*DEXP(ED*TX) 

ELSE 
END  IF 

IF(J5.EQ.27)  THEN 
XI  =  AO*DEXP(-(EA+DELl)*TX) 

X2  =  BO*DEXP(EB*TX) 

X3  =  D0*DEXP((ED+DEL4)*TX) 

ELSE 

ENDIF 

IF(J5.EQ.28)  THEN 
XI  =  AO*DEXP(-EA*TX) 

X2  =  B0*DEXP((EB+DEL2)*TX) 

X3  =  D0*DEXP((ED+DEL4)*TX) 

ELSE 

ENDIF 


setting  initial  conditions  at  reactor  entrance 


Y0(1)  =  1.0D00 


EPISODE  error  control 

IERROR  =  1  Absolute  error  control,  Ymax(i)  =  1 
IERROR  =  2  Error  relative  to  Y(i)  is  controlled. 

If  Y(i)  =  0,  a  divide  error  occurs 
IERROR  =  3  Error  relative  to  the  largest  value  of 
Y(i)  is  controlled. 


IERROR  =  1 
TO  =  0.0D00 
MF  =  22 
INDEX  =  1 
HO  =  SAVEH 

CALL  DRIVE(N,  TO,  HO,  Y0,  TOUT,  EPS,  IERROR,  MF,  INDEX) 
CONX(J5)  =  (ONE  -  Y0(1))*100.0D00 
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IF(J5.EQ.  1)  THEN 
X(J)  =  CONX(l) 

\VRITE(*,129)  X(J),  CONV(J) 

WRITE(9,129)  X(J),  CONV(J) 

SUM  =  SUM  +  (X(J)-CONV(J))*(X(J)-CONV(J)) 

ELSE 

ENDIF 

400  CONTINUE 

PA  =  (CONX(2)  -  CONX(l))/DELA 
PB  =  (CONX(4)  -  CONX(l))/DELB 
PD  =  (CONX(6)  -  CONX(l))/DELD 
PI  =  (CONX(8)  -  CONX(I))/DELl 
P2  =  (CONX(IO)  -  CONX(l))/DEL2 
P4  =  (CONX(12)  -  CONX(l))/DEL4 

PAA  =  (CONX(2)  -  TWO*CONX(l)  +  CONX(3))/(DELA*DELA) 

PBB  =  (CONX(4)  -  TWO*CONX(l)  +  CONX(5))/(DELB*DELB) 

PDD  =  (CONX(6)  -  TWO*CONX(l)  +  CONX(7))/(DELD*DELD) 

PH  =  (CONX(8)  -  TWO*CONX(l)  +  C0NX(9))/(DEL1*DEL1) 

P22  =  (CONX(IO)  -  TWO*CONX(l)  +  CONX(ll))/(DEL2*DEL2) 

P44  =  (CONX(12)  -  TWO*CONX(l)  +  CONX(13))/(DEL4*DEL4) 

PAB  =  (CONX(l)  -  CONX(2)  -  CONX(4)  +  CONX(14))/(DELA*DELB) 
PAD  =  (CONX(l)  -  CONX(2)  -  CONX(6)  +  CONX(15))/(DELA*DELD) 
PA1  =  (CONX(l)  -  CONX(2)  -  CONX(8)  +  CONX(16))/(DELA*DELl) 
PA2  =  (CONX(l)  -  CONX(2)  -  CONX(IO)  +  CONX( 1 7))/(DELA*DEL2) 
PA4  =  (CONX(l)  -  CONX(2)  -  CONX(12)  +  CONX(I8))/(DELA*DEL4) 
PBD  =  (CONX(l)  -  CONX(4)  -  CONX(6)  +  CONX( 1 9))/(DELB*DELD) 
PB1  =  (CONX(l)  -  CONX(4)  -  CONX(8)  +  CONX(20))/(DELB*DEL1) 
PB2  =  (CONX(l)  -  CONX(4)  -  CONX(IO)  +  CONX(21))/(DELB*DEL2) 
PB4  =  (CONX(l)  -  CONX(4)  -  CONX(12)  +  CONX(22))/(DELB*DEL4) 
PD1  =  (CONX(l)  -  CONX(6)  -  CONX(8)  +  C0NX(23))/(DELD*DEL1) 
PD2  =  (CONX(l)  -  CONX(6)  -  CONX(IO)  +  CONX(24))/(DELD*DEL2) 
PD4  =  (CONX(I)  -  CONX(6)  -  CONX(12)  +  CONX(25))/(DELD*DEL4) 
P12  =  (CONX(l)  -  CONX(8)  -  CONX(IO)  +  CONX(26))/(DELl*DEL2) 
P14  =  (CONX(l)  -  CONX(8)  -  CONX(12)  +  C0NX(27))/(DEL1*DEL4) 
P24  =  (CONX(l)  -  CONX(IO)  -  CONX(12)  +  CONX(28))/(DEL2*DEL4) 

FI  =  FI  +  TWO*CONV(J)*PA  -  TWO*CONX(l)*PA 
F2  =  F2  +  TWO*CONV(J)*PB  -  TWO*CONX(l)*PB 
F3  =  F3  +  TWO*CONV(J)*PD  -  TWO*CONX(l)*PD 
F4  =  F4  +  TWO*CONV(J)*Pl  -  TWO*CONX(l)*PI 
F5  =  F5  +  TWO*CONV(J)*P2  -  TWO*CONX(l)*P2 
F6  =  F6  +  TWO*CONV(J)*P4  -  TWO*CONX(l)*P4 

F1A  =  F1A  +  TWO*(CONV(J)*PAA  -  CONX(l)*PAA  -PA*PA) 

FIB  =  FIB  +  TWO*(CONV(J)*PAB  -  CONX(l)*PAB  -PA*PB) 

FID  =  FID  +  TWO*(CONV(J)*PAD  -  CONX(l)*PAD  -PA*PD) 

FI  1  =  FI  1  +  TWO*(CONV(J)*PAl  -  CONX(l)*PAl  -PA*P1) 

F12  =  F12  +  TWO*(CONV(J)*PA2  -  CONX(I)*PA2  -PA*P2) 

F14  =  F14  +  TWO*(CONV(J)*PA4  -  CONX(l)*PA4  -PA*P4) 

F2A  =  F2A  +  TWO*(CONV(J)*PAB  -  CONX(l)*PAB  -PB*PA) 

F2B  =  F2B  +  TWO*(CONV(J)*PBB  -  CONX(l)*PBB  -PB*PB) 

F2D  =  F2D  +  TWO*(CONV(J)*PBD  -  CONX(l)*PBD  -PB*PD) 

F2I  =  F21  +TWO*(CONV(J)*PBl  -C0NX(1)*PB1  -PB*P1) 
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F22  =  F22  +  TW0*(C0NV(J)*PB2  -  C0NX(1)*PB2  -PB*P2) 
F24  =  F24  +  TW0*(C0NV(J)*PB4  -  C0NX(1)*PB4  -PB*P4) 

F3A  =  F3A  +  TWO*(CONV(J)*PAD  -  C0NX(1)*PAD  -PD*PA) 
F3B  =  F3B  +  TWO*(CONV(J)*PBD  -  C0NX(1)*PBD  -PD*PB) 
F3D  =  F3D  +  TWO*(CONV(J)*PDD  -  C0NX(1)*PDD  -PD*PD) 
F3 1  =  F3 1  +  TWO*(CONV(J)*PD  1  -  C0NX(1)*PD1  -PD*P1) 
F32  =  F32  +  TW0*(C0NV(J)*PD2  -  C0NX(1)*PD2  -PD*P2) 
F34  =  F34  +  TW0*(C0NV(J)*PD4  -  C0NX(1)*PD4  -PD*P4) 

F4A  =  F4A  +  TW0*(C0NV(J)*PA1  -  C0NX(1)*PA1  -P1*PA) 
F4B  =  F4B  +  TW0*(C0NV(J)*PB1  -  C0NX(1)*PB1  -Pl+PB) 
F4D  =  F4D  +  TWO*(CONV(J)*PD  1  -  C0NX(1)*PD1  -P1*PD) 
F41  =F41  +  TW0*(C0NV(J)*P11  -C0NX(1)*P11  -P1*P1) 

F42  =  F42  +  TW0*(C0NV(J)*P12  -  C0NX(1)*P12  -P1*P2) 

F44  =  F44  +  TW0*(C0NV(J)*P14  -  C0NX(1)*P14  -P1*P4) 

F5A  =  F5A  +  TW0*(C0NV(J)*PA2  -  C0NX(1)*PA2  -P2*PA) 
F5B  =  F5B  +  TW0*(C0NV(J)*PB2  -  C0NX(1)*PB2  -P2*PB) 
F5D  =  F5D  +  TW0*(C0NV(J)*PD2  -  C0NX(1)*PD2  -P2*PD) 
F51  =  F51  +  TW0*(C0NV(J)*P12  -  C0NX(1)*P12  -P2*P1) 

F52  =  F52  +  TWO*(CONV(J)*P22  -  CONX(l)*P22  -P2*P2) 

F54  =  F54  +  TWO*(CONV(J)*P24  -  CONX(l)*P24  -P2*P4) 

F6A  =  F6A  +  TW0*(C0NV(J)*PA4  -  C0NX(1)*PA4  -P4*PA) 
F6B  =  F6B  +  TW0*(C0NV(J)*PB4  -  C0NX(1)*PB4  -P4*PB) 
F6D  =  F6D  +  TW0*(C0NV(J)*PD4  -  C0NX(1)*PD4  -P4*PD) 
F61  =  F61  +  TW0*(C0NV(J)*P14  -  C0NX(1)*P14  -P4*P1) 

F62  =  F62  +  TWO*(CONV(J)*P24  -  CONX(l)*P24  -P4*P2) 

F64  =  F64  +  TWO*(CONV(J)*P44  -  CONX(l)*P44  -P4*P4) 

500  CONTINUE 
100  CONTINUE 
Al(l,7)  =  -FI 
A  1(2 ,7)  =  -F2 
Al(3,7)  =  -F3 
Al(4,7)  =  -F4 
Al(5,7)  =  -F5 
Al(6,7)  =  -F6 
A1(1,1)  =  F1A 
Al(l,2)  =  FIB 
Al(l,3)  =  FID 
Al(l,4)  =  FI  1 
Al(l,5)  =  F12 
Al(l,6)  =  F14 
A  1(2,1)  =  F2A 
A  1(2,2)  =  F2B 
Al(2,3)  =  F2D 
A  1(2,4)  =  F21 
Al(2,5)  =  F22 
Al(2,6)  =  F24 
Al(3,l)  =  F3A 
Al(3,2)  =  F3B 
Al(3,3)  =  F3D 
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Al(3,4)  =  F31 
Al(3,5)  =  F32 
Al(3,6)  =  F34 
Al(4,l)  =  F4A 
Al(4,2)  =  F4B 
A  1(4,3)  =  F4D 
A  1(4,4)  =  F41 
Al(4,5)  =  F42 
A  1(4,6)  =  F44 
Al(5,l)  =  F5A 
Al(5,2)  =  F5B 
Al(5,3)  =  F5D 
A  1(5,4)  =  F51 
Al(5,5)  =  F52 
Al(5,6)  =  F54 
A1(6,1)  =  F6A 
Al(6,2)  =  F6B 
A  1(6,3)  =  F6D 
Al(6,4)  =  F61 
Al(6,5)  =  F62 
Al(6,6)  =  F64 
206  format(7dl0.4) 

CALL  GAUSS(KORD,  KORD1,  Al,  P5) 

AOLD  =  AO 
BOLD  =  BO 
DOLD  =  DO 
OLD 1  =  E  A 
OLD2  =  EB 
OLD4  =  ED 
AO  =  AO  +  P5(l) 

BO  =  BO  +  P5(2) 

DO  =  DO  +  P5(3) 

EA  =  EA  +  P5(4) 

EB  =  EB  +  P5(5) 

ED  =  ED  +  P5(6) 

EAO  =  ((AOLD  -  AO)/AO)*((AOLD  -  AO)/ AO) 
EBO  =  ((BOLD  -  BO)/BO)*((BOLD  -  BO)/BO) 
EDO  =  ((DOLD  -  DO)/DO)*((DOLD  -  D0)/D0) 
EEA  =  ((OLD  1  -  EA)/EA)*((OLD  1  -  EA)/EA) 
EEB  =  ((OLD2  -  EB)/EB)*((OLD2  -  EB)/EB) 
EED  =  ((OLD4  -  ED)/ED)*((OLD4  -  ED)/ED) 
IF(AO.gt.2,ODOO*AOLD)  AO  =  2.0D00*AOLD 
IF(A0.1t.0.5D00*AOLD)  AO  =  0.5D00*AOLD 
IF(B0.gt.2.0D00*BOLD)  BO  =  2,ODOO*BOLD 
IF(B0.1t.0.5D00*BOLD)  BO  =  0.5D00*BOLD 
IF(DO.gt.2,ODOO*DOLD)  DO  =  2.0D00*DOLD 
IF(D0.1t.0.5D00*DOLD)  DO  =  0.5D00*DOLD 
IF(EA.gt.  1.2DOO*OLD1)  EA  =  1.2DOO*OLD1 
IF(EA.lt.0.8D00*OLDl)  EA  =  0.8D00*OLDl 
IF(EB.gt.  1.2DOO*OLD2)  EB  =  1.2DOO*OLD2 
IF(EB.lt.O,8DOO*OLD2)  EB  =  0.8D00*OLD2 
IF(ED.gt.l.2DOO*OLD4)  ED  =  1.2DOO*OLD4 
IF(ED.lt.O,8DOO*OLD4)  ED  =  0.8D00*OLD4 
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SUME  =  dsqrt(EA0+EB0+ED0+EEA+EEB+EED)/5.0D00 
ERROR  =  DSQRT(SUM)/DBLE(NDAT-1) 

WRITE/*,  551)  ERROR 
WRITE(9,  551)  ERROR 

WRITER,  547)  AO,  BO,  DO,  EA,  EB,  ED,  SUME 
WRITE(9,  547)  AO,  BO,  DO,  EA,  EB,  ED,  SUME 
SUM  =  0.0D00 

IF(SUME.LT.1.0D-04)  GOTO  699 
600  CONTINUE 
699  CONTINUE 
800  CONTINUE 

547  FORMAT(/, 2X,' AO  =  ’,D12.6,2X,'B0  =  ',D12.6,2x,'D0  =  ',D12. 
1  2X,'EA  =  ',D12.6,2X,'EB  =  ',D12.6,2x,'ED  =  ',D12.6,/, 

1  2x, 'ERROR  =  ',dl2.6,/) 

551  FORMAT(/2X,'XNORM  =  ’,F8.5) 

119  FORMAT(2x, ’Conversion: ’,fl0.4) 

129  FORMAT(2x,'Pred.  Conv:  ’,fl0.4,2x,'Exp.  Conv:  ’,fl0.4) 

999  CONTINUE 
close(9) 

STOP 

END 

SUBROUTINE  GAUSS(KORD,KORDl,A,P5) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  A(6,7),  P5(7) 

KORD  =  6 
KORD 1  =  KORD  +  1 
KORD2  =  KORD  -  1 
DO  100  N1  =  l,KORD 
N3  =N1  +  1 
DO  101  N2  =  N3,KORD 

IF(DABS(A(N2,N1)).GT.DABS(A(N1,N1)))  GOTO  105 
GOTO  101 

105  DO  106  I  =  1.KORD1 
SAVE1  =  A(N2,I) 

SAVE2  =  A(N1,I) 

A(N1,I)  =  SAVE1 

106  A(N2,I)  =  SAVE2 
101  CONTINUE 

DO  115  N  =  N1, KORD 
ASAVE  =  A(N,N1) 

IF(DABS(ASAVE).LE.1.0D-50)  GOTO  115 
DO  1 16  M  =  Nl,KORDl 
A(N,M)  =  A(N,M)/ASAVE 
116  CONTINUE 
115  CONTINUE 

IF(Nl.EQ.KORD)  GOTO  100 
NEQ  =  N1  +  1 
DO  120  N  =  NEQ, KORD 
IF(D ABS(A(N,N  1  )).LE.  1 .0D-50)  GOTO  120 
DO  121  N2  =  Nl,KORDl 
IF(DABS(A(N1,N)).LE.  1.0D-50)  GOTO  121 
A(N,N2)  =  A(N1,N2)  -  A(N,N2) 

121  CONTINUE 


APPENDIX 


38 


n  n  o  o 


120  CONTINUE 
100  CONTINUE 

P5(KORD)  =  A(KORD,KORDl) 

DO  200  N=  1.KORD2 
RHS  =  0.0D00 
DO  210  I  =  1,N 

210  RHS  =  P5  (KORD 1  -I)  *  A(KORD  -N.KORD  1  -I)  +  RHS 

200  P5(KORD-N)  =  A(KORD-N,KORD  1 )  -  RHS 
RETURN 
END 

SUBROUTINE  DIFFUN  (N,  T,  Y,  YDOT) 

IMPLICIT  REAL*8  (A-H,  O-Z) 

DIMENSION  Y(N),  YDOT(N) 

COMMON  /SBCOM1/XKM1,  SA,  SS,  CA0,  RHOB 
COMMON  /SBCOM3/  XI,  X2,  X3 
COMMON  /SBCOM2/  WASHT,  CELLR,  TOUT,  VEL 
common  /sbcom3/  j,  j5 

DATA  PT5  /0.5D00/,  ONE  /1.0D00/,  TWO  /2.0D00/,  FOUR/4.0D00/ 


c  The  surface  concentration  of  HCN,  XAS,  is  determined  using 
c  Newton's  method  of  itteration,  where 
c 

c  0  =  XKM1*SA*(Y(1)  -XAS)-  Rate  =  f(XAS) 
c 

XAS  =  Y(l) 

15  FORMAT(2X,'XAS  CONVERGENCE  FAILURE,  XAS  =  ',D12.6,/, 
1  2X, 'ERROR  =  ',D  12.6) 


SS  surface  area,  cm**2  catalyst/cm* *3  reactor  void 
SA  surface  area,  cm* *2  catalyst/cm**3  catalyst 


A1  =  -XKM1*SA*CA0*CA0*X3 
B1  =  X3*CA0*CA0*XKM1*SA*Y(1)  -  XKM1*SA*CA0  - 
1  XKM1*SA*CAO*CAO*X2*(ONE  -  Y(l))  -  X1*CA0 
Cl  =XKM1*SA*CA0*Y(1)  +  XKM1*SA*CA0*CA0*X2*(Y(1)  -  Y(1)*Y(1)) 
XAS  =  (-B1  -  DSQRT(B1*B1  -  FOUR*Al*Cl))/(TWO*Al) 

IF(XAS.GT.ONE)  WRITE(*,15)  XAS 

YDOT(l)  «  -  (0NE/VEL)*XKM1*SS*(Y(1)  -  XAS) 

RETURN 

END 

SUBROUTINE  PEDERV(N,T,Y,PD,NO) 

IMPLICIT  REAL*8  (A-H,  O-Z) 

DIMENSION  Y(N0),  PD(N0,N0) 

RETURN 

END 


APPENDIX 


39 


nnnn 


File  M0N01.DAT  provides  input  parameters  to  run 
Fit  parameter  routine,  MONOl.F 


52 

.46429200E+04 

.46250500E+08 

.33751 100E+08 

.60752800E+04 

.44557700E+04 

.26948300E+04 

.10000000E+01 

.10000000E-04 

.50000000E+01 

.10000000E-05 

.57315000E+03 

.13061000E+01 

.10370000E+00 

.35000000E-02 

.1 1500000E-01 

.77500000E+01 

.50000000E+04 

.50000000E+04 

.50000000E+04 

.10000000E+02 

.10000000E+02 

.10000000E+02 


NUMBER  OF  DATA  POINTS 
v.  AO 
BO 
DO 
EA/R 
EB/R 
ED/R 

CATALYST  BED  DENSITY  (G/CM**3) 
INITIAL  STEP  SIZE  (CM) 

Catalyst  Length,  cm  (TOUT) 

ERROR  TOLERANCE  (EPS) 

Reference  Temperature,  K  (TREF) 

PRESSURE  (ATM)  (PRESS) 

CELL  RADIUS  (CM) 

WASHCOAT  THICKNESS 

SUBSTRATE  THICKNESS 

MASS  TRNASFER  COEFFICIENT,  CM**2/S 

DELA 

DELB 

DELD 

DELI 

DEL2 

DEL3 
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DISTRIBUTION  LIST  FOR  ERDEC-CR-219 


NAMES  COPIES 

CDR  USA  CBDCOM 

AMSCB  CG  E4435  UPCHURCH  1 
AMSCB  CII  E3330  LIBRARY  2 
AMSCB  CIH  E5183  SMART  1 
5232  FLEMING  ROAD 


ABERDEEN  PROVING  GROUND  MD  21010-5423 


DIR  USA  ERDEC 

SCBRD  ASI  E5101  FAMINI  1 
SCBRD  RT  E3330  TECH  REL  OFC  1 
SCBRD  RT  E3330  J  WILLIAMS  1 
SCBRD  RTE  E3330  D  REED  25 
5232  FLEMING  ROAD 


ABERDEEN  PROVING  GROUND  MD  21010-5423 
DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS 
ACTIVITY 

ATTN  AMXSY  CB  W  HEAPS  1 

392  HOPKINS  ROAD 

ABERDEEN  PROVING  GROUND  MD  21005-5071 
COMMANDANT 

US  ARMY  INFANTRY  SCHOOL 
ATTN  ATSH  WC 
FORT  BENNING  GA  31905-5400 


NAMES  COPIES 


COMMANDER 

US  ARMY  ARMOR  CENTER 

MOUNTED  WARFIGHTING  BATTLES PACE  BATTLE 
LABORATORY 

ATTN  ATZK  MW  1 

FORT  KNOX  KY  40121-5000 

COMMANDER 

US  ARMY  READINESS  GROUP 

ATTN  AFKB  RG  FS  CS  1 

FORT  SILL  OK  73503-6700 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
8725  JOHN  J  KINGMAN  ROAD  12 

SUITE  0944 

FT  BELVOIR  VA  22060-6218 
OUSD/DTSA/TD 

ATTN  PATRICIA  SLYGH  1 

400  ARMY  NAVY  DRIVE  ROOM  305 
ARLINGTON  VA  22202 


1 


